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We study a zero-dimensional version of the fluctuating nonlinear hydrodynamics (FNH) of super- 
cooled liquids originally investigated by Das and Mazenko (DM) [Phys. Rev. A 34, 2265 (1986)]. 
The time-dependent density-like and momentum-like variables are introduced with no spatial de- 
grees of freedom in this toy model. The structure of nonlinearities takes the similar form to the 
original FNH, which allows one to study in a simpler setting the issues raised recently regarding the 
field theoretical approaches to glass forming liquids. We study the effects of density nonlinearities on 
the time evolution of correlation and response functions by developing field theoretic formulations 
in two different ways: first by following the original prescription of DM and then by constructing a 
dynamical action which possesses a linear time reversal symmetry as proposed recently. We show 
explicitly that, at the one- loop order of the perturbation theory, the DM- type field theory does not 
support a sharp ergodic-nonergodic transition, while the other admits one. The simple nature of 
the toy model in the DM formulation allows us to develop numerical solutions to a complete set of 
coupled dynamical equations for the correlation and response functions at the one-loop order. 
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I. INTRODUCTION 



The slow dynamics of supercooled liquids near the 
glass transition has been under intense theoretical and 
experimental investigation for many years. Among many 
theoretical attempts to understand the slowing down of 
supercooled liquids, the mode coupling theory (MCT) 
P) S H, 01 stands out as one of the most successful ones. 
It explains, for example, an elaborate sequence of time 
relaxation processes with characteristic exponents which 
are consistent with experimental findings. In its initial 
form d, 1^, referred to as the standard MCT, it pre- 
dicts a sharp ergodic-to-nonergodic (ENE) transition at a 
critical temperature or density with the nonergodic phase 
characterized by the density autocorrelation function ap- 
proaching a nonzero value called the nonergodicity pa- 
rameter in the long time limit. Many experiments and 
numerical simulations, however, show that this feature is 
not realized in real supercooled liquids in finite dimen- 
sions, and that the ergodicity is restored at finite tem- 
perature. 

There have been many atternpts to put the MCT into 
a field theoretic framework d, M, M, [ill, [S M, Q M, 
m, 113 , since it has many advantages including the pos- 
sibility of a systematic improvement. Das and Mazenko 
(DM) [8] studied the nonlinear feedback mechanism of 
density fluctuations in supercooled liquids by formulat- 
ing a field theoretic renormalized perturbation theory of 
the fluctuating nonlinear hydrodynamics (FNH) of com- 
pressible fluids. They find that the sharp transition is 
cutoff and the system remains ergodic at all tempera- 
tures or densities. Recently, however, the validity of the 
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DM results, especially of those on the explanation of the 
cutoff mechanism, was questioned in Refs. [Til, [13, fisj . 
Some of these works are based on a field theory devel- 
oped in Ref. TT| where the dynamical action is invariant 
under a set of linear time-reversal transformations. This 
formulation allows one to have a full set of fluctuation- 
dissipation relations (FDR) relating linearly correlation 
functions to response functions. The fleld theory of DM 
has only a limited number of linear FDR and some re- 
lations hold only in the hydrodynamic limit. This is 
one of the points on which the conclusion by DM on 
the absence of the sharp ENE transition was questioned. 
The field theory with linear FDR was later improved [l3| 
for the case of interacting Browni an p articles satisfying 
the Dean-Kawasaki equation [l^ |20| |. where the stan- 
dard MCT result was recovered at the one-loop order of 
perturbation theory. This improved method was then 
applied to the FNH [l7| with results indicating a sharp 
ENE transition at the one-loop order with the noner- 
godicity parameter satisfying the standard MCT result. 
In response to these developments, DM reexamined their 
work and showed [2]| in a nonperturbative analysis with- 
out resorting to the hydrodynamic limit that the sharp 
ENE transition is not present in the FNH after all. This 
conclusion is also supported by the recent direct numer- 
ical integration of the generalized Langevin equations of 
the FNH [22]. 

It is somewhat puzzling to have completely different re- 
sults from the two field theoretic approaches of the same 
model. In this respect, it might be desirable to have 
a simpler setting in which one can compare these two 
field theoretical approaches and study where the differ- 
ence originates. The field theoretical treatment of FNH 
is complicated by many factors including the presence of 
many dynamical variables. This is especially the case for 
the field theory with linear FDR. In this paper, we intro- 
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duce a simple toy model of FNH, which can shed some 
light on the issues described above concerning the field 
theoretic approaches to the FNH. In this toy model, there 
is no spatial dependence in the dynamical field variables 
which consist simply of a density-like variable and a single 
component momentum-like variable. We develop the two 
different field theories of the toy model, namely the orig- 
inal DM-type field theory and the one with linear FDR. 
We show explicitly that, at the one- loop order of pertur- 
bation expansion, a sharp ENE-type transition does not 
occur in the DM-type field theory. On the other hand, 
the field theory with linear FDR shows a sharp transi- 
tion at the one-loop order. By comparing the two field 
theories, we find that the major difference between the 
two formulations lies in the way of treating the density 
nonlinearities present in FNH within the renormalized 
perturbation theory. In particular, the field theory with 
linear FDR results in a dynamical action which contains 
non-polynomial functions of field variables in contrast to 
the DM field theory. This implies that, when the renor- 
malized perturbation theory is performed at a given order 
of the loop expansion, the two field theories end up with 
treating the density nonlinearities in a different way, since 
the field theory with linear FDR involves truncating the 
non-polynomial functions. 

Although there is only a limited number of linear FDR, 
the DM field theory at a given order of the loop expansion 
can be regarded as a well-defined self-consistent theory 
among the correlation and response functions satisfying a 
set of coupled self-consistent equations. In this paper, we 
construct the set of coupled equations at the one-loop or- 
der for the DM field theoretic approach to the toy model 
and study them numerically. Since the simple nature of 
the toy model reduces the number of independent corre- 
lation and response functions, we were able to solve these 
equations numerically. 

In the next section, we introduce our toy model and 
construct the field theories following the DM prescrip- 
tion and the method involving linear FDR, respectively. 
In Sec. mil we study the time evolution of the correla- 
tion functions using the Schwinger-Dyson equations for 
both field theoretic formulations. We then analyze the 
time evolution equations for the correlation function cor- 
responding to the density autocorrelation function of the 
FNH at the one-loop order for the possible existence of 
a sharp ENE-type transition. In Sec. IIVI we present a 
set of coupled equations for the correlation and response 
functions in the DM field theory at the one-loop order 
and their numerical solutions. In the final section, we 
summarize our results with discussion. 



II. MODEL 



component momentum-like variable b(t) without any spa- 
tial degrees of freedom. In order to construct the equa- 
tions motion for these variables, we introduce the effec- 
tive free energy F where the equilibrium distribution for 
the system at temperature T is given by exp{—F/T). 
The free energy can be written as -F = + Fjj where 
FK[a.,b] is the kinetic energy and the potential energy 
part -^[/[a] is assumed to depend only on a. We take the 
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usual form for the kinetic part, that is Fk — The 
nonlincarity in the form of 1/a plays an important role 
in the following discussion. 

Kim and Kawasaki [23j introduced a similar zero- 
dimensional toy model of the FNH to the present one 
by incorporating the multicomponent density-like and 
momentum-like variables. However, in addition to having 
multicomponent fields, their model differs from ours in a 
fundamental way. Their free energy is quadratic both in 
the density-like and the momentum-like variables with- 
out the 1/a nonlincarity which is present in the FNH of 
compressible fluids. In the present toy model, therefore, 
the actual form of the equations of motion will be dif- 
ferent from those in Ref. [l^, but the derivation of the 
equations from the free energy can be performed in the 
same way. In our model, the equation of motion for the 
variable a{t) takes the form of a zero-dimensional version 
of the continuity equation, namely 



a{t) + Jb{t) = 



(1) 



for some constant J. This can be regarded as the re- 
versible dynamics for a{t) which can be derived from 



ait) 



dF dQab 

ah— J - 



db 



db 



(2) 



with Qab ^ —Ja playing the role of the Poisson bracket. 
The equation of motion for b(t) has the dissipative part 
described by the coefficient F in addition to the reversible 
part as follows: 



bit) 



^ dF dQba 
>ba— J - 



da 



da 



dF 



(3) 



where Qba = ~Qab and the Gaussian white noise 9 has 
the variance (6'(t)6l(t')) = 2rTSit - t'). In the present 
work, we take the simple quadratic form for the potential 
energy part Fu of the effective free energy F as 



Fu[a] = ^iSa)' 



(4) 



with the fluctuation Sa — a — ao and the average value 
aQ. We thus have the equation of motion for b{t) as 



bit) + J{—^- JAaiSa) +TJ + t(- 
2a I \a 



(5) 



Our model is a zero-dimensional version of the FNH 
of compressible fluids developed by Das and Mazenko 
Q. We introduce as our dynamical variables a 
time-dependent density-like variable a(i) and a single- 



We can easily verify that the equilibrium stationary 
distribution corresponding to the above equations for a 
and b is proportional to exp(— F/T). The corresponding 
Fokker-Planck equation for the probability distribution 
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P{a, b, t) is given by dtP = CP where the Fokker-Planck 
operator is given by £ = £i + £2 with 



db \ db a 



(6) 



JAa{Sa) + TJ (7) 



It is straightforward to show that P exp(— F/T) satis- 
fies CP = 0. 

One can develop a field theory from the above 
Langevin equations by using the standard Martin-Siggia- 
Rose (MSR) formahsm In the MSR procedure, the 
hatted fields a and b are introduced to enforce the equa- 
tions of motion for a and b respectively. We present 
in the following subsections two different field theoret- 
ical approaches to this model and compare the outcomes 
of both approaches concerning the existence of an ENE 
transition. The first one is the original Das and Mazcnko 
approach 8], where the 1/a nonlinerarities in the model 
are taken care of in a simple way by the introduction of 
a single additional auxiliary field. The second approach 
due to Refs. [Ill Hi [13 incorporates the linear time- 
reversal symmetry into the dynamical action resulting in 
a set of linear FDR. In order to do that more auxiliary 
fields will have to be introduced. 



A. Das-Mazenko approach 

In the DM approach, an auxiliary velocity-like field 
c(t) is introduced such that the condition b{t) — a[t)c[t) 
is enforced through a delta function 



1 = / Vc{t) 5{a{t)c{t) - b{t)) 

Vc{t) / Vc{t) exp [-ic{t) (6(t) - a{t)c{t))] .(8) 



The first equality holds up to a Jacobian. This Jacobian 
was shown in Ref. |25j to have no effect on the corre- 
lation and response functions and will be neglected in 
the following analysis. Using this identity, we obtain the 
generating functional Z as a functional integral over the 
fields ipi{t) = 5a{t),b(t),c{t) and ^i{t) = a{t) ,b{t) , c{t) . 

We can write Zdm = / ^V'j^^V'j exp(-S'DM[V', V']), 
where 



dt [TTb^{t) + ia{t) {d{t) + Jb{t)} 



+ib{t){bit) - JAaoSa{t) + Tc{t) + TJ 

+ ^{c{t)f + ^6a{t){c{t)f JA{6a{t)f} 



+ic{t){b{t) - aoc{t) ~ Sa{t)c{t)}]. 



(9) 



We will use \E'(i) to represent any one of the six vari- 
ables {tpi,ipi} in our model, and denote the two-point 



correlation function between arbitrary two variables ^(t) 
and -^'{t') by 



(10) 



(For the subscripts of G, we will use a instead of Sa for 
simplicity.) Note that among the correlation functions 
those between two hatted variables vanish due to causal- 
ity, that is = 0. It follows that iTJb{t) term in the 
action Eq. ^ has no effect on the correlation functions 
and will be neglected in the following. The causality also 
requires that G^,^ \t) — Q for t < 0. 

We can easily establish some nonperturbative relations 
among the correlation functions which will be useful in 
later discussion. If we use ^(i') — ^Ji{t') = Sa{t'),b{t') 
and c{t') in the identity 



= JUv^.Vi;, ^ [^{t') exp(-5DM)] , 



we obtain 



Ga.^it) + JGi,^it)=0. 



(11) 



(12) 



On the other hand, if *(t') 
are used, we have 



TP{t') ^ a(t'),6(t') and c(t') 



Gadt) + JG,.{t) = -^^r^^{t). 



(13) 



Note that Eqs. (fT2|) and ([T3|) are direct consequences 
of the zero-dimensional version of the mass conservation 
law given by Eq. ([T]). In the DM approach, only a lim- 
ited number of FDR exist that relate linearly the cor- 
relation functions to response functions. Assuming the 
time reversal properties of the fields as a(—t) — a{t), 
b{—t) = —b{t) and c{—t) — ~c{t), we can derive the 
FDR for = a,b and c as 



--e(i)G^e(i), 



(14) 



where 0(t) = 1 for t > and vanishes for t < 0. The de- 
tailed derivation of the FDR closely follows the one given 
in Ref. [8|. Since ipi is a real field, we can show that the 
correlation function between unhatted and hatted vari- 
ables is a pure imaginary number, that is 



G 



(15) 



B. Field theory with linear FDR 

We apply the field theoretical approach developed in 
Refs. [nl, [ij) [13 to our toy model given by the dynamic 
equations (U) and Among these methods, we will fol- 
low closely the improved procedure described in Ref. jl7| . 
In order to do that, we introduce two auxiliary fields r]{t) 



and 9{t) defined by 



ao ao f-j' 



db 
f,j{6a,b) 

.,c ^ b'^ v^/ ( ^a 



t_ 



fc=0 



ao 



(17) 



Note that these are non-polynomial functions of the main 
dynamical variables 5a and b. Introducing the hatted 
counterparts f}{t) and 9{t) to enforce these definitions 
for the new variables, we can construct the generating 
functional as Zfdr = J lli'^4>i'^(t>i exp(-S'FDR[0, 0]), 
where <p,{t) = da{t),b{t),r]{t),0{t) and 4)^{t) = 
a{t),b{t),'ri{t),6{t). Similarly to the DM case, we will 
use $(<) to represent one of the eight variables {(j)i,(t>i}. 
The dynamical action in this case can be written as a 
sum of the Gaussian and the nonlinear parts as Sfdr — 
•^FDR + "^FDR' where 



5, 



FDR = / dt [TTb\t) + ia{t){a{t) + Jb{t) + Jaor^{t)} 

+ib{t){b{t) - JAaoSa{t) - Jaoe{t) + rc{t) + —b{t)} 

ao 



-tfi{t)'n{t) + id{t)9{t)] 



(18) 



and 



4dr = / dt [ia{t){JSa{t)T]{t) + —Sa(t)b{t)} 



J_ 

ao 



-ib{t){JSa{t)9{t) + JA{Sa{t)f} 
-ifiit)USa,b)~i9it)fgi6a,b)]. (19) 

Note that we have used the identity 

1 



aoc(t) + Sa{t)T]{t) + —Sa{t)b{t) = 0, 
ao 



(20) 



which follows directly from the definition of r]{t), Eg. p^ . 

The above actions ■S'ppf^ and ■S'ppj^ are separately invari- 
ant under the time reversal transformations given by 

Sa{-t) = da{t), b{-t) = ~b{t), 
Tji-t) = ~rj{t), 9{-t) = 9{t), 



a{-t) 
bi-t) : 
fji-t) 
9{-t) 



--a{t)'^m-^pa{t), 

^ -m + ^m, 
m-^a{t). 



(21) 



Applying the time reversal invariance on the various cor- 
relation functions by following the procedures described 





(16) FIG. 1: The one-loop diagrams contributing to the self-energy 



in Refs. [ll|, IM ll7|, we obtain a set of FDR. Here we 
only list those which are relevant to the discussion in the 
next section. We have for (p — a,b, rj and 9 



G^a{t) — 


-^e(t)[AG^,(0-f 


' Gt^g{t)\, 


G^iit) = 


-^mi-GMt)- 

T ao 


f G^,(i)], 


G^fiii) = 






G^eit) = 







(22) 



III. 



RENORMALIZED PERTURBATION 
THEORY: ONE-LOOP ORDER 



In this section, we develop self-consistent renormal- 
ized perturbation theories for the two field theoretic ap- 
proaches introduced in the previous section. We then 
focus on the time evolution of two-point correlation func- 
tions using the Schwinger-Dyson (SD) equation. In par- 
ticular, we study the t ^ oo limit of Gaa (t) , which corre- 
sponds to the density auto-correlation function in FNH, 
to explore the possibility of an ENE transition. The for- 
mal development of the self-consistent perturbation the- 
ory can be found in Refs. [Illli, [23|. The SD equation 
defines the self-energy Y. through its relation to the prop- 
agator (the two-point correlation function) G. It is given 
symbolically by 



(23) 



where the subscript refers to the bare quantity ob- 
tained by keeping only the Gaussian terms in the ac- 
tion. The self-energy is obtained by differentiating the 
so-called two-particle irreducible vertex function r2Pi[G] 
with respect to the propagator G. At the one- loop order 
of the loop expansion of r2Pi, there are only two kinds 
of diagrams for the self-energy which are relevant to the 
two field theoretical approaches studied in the previous 
section. These are shown in Fig. [TJ 
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A. Das-Mazenko approach 

The SD equation for the DM field theory between ar- 
bitrary two fields 5* and is given by 

(5>i(*'(5(t — t') 

= E / dt"{[G^']^^„ it - no.,,,.,, it" - 1') 

-t 

it-t")G^,,^>it" ~t')}. (24) 

We note that the causality requirement for the self-energy 
reads T,^,^\t) = for i < and E^.^^. (t) = for all 
t. Similarly to the propagators, the self-energy between 
unhatted and hatted variables is pure imaginary. Since 
there is no nonlinear term containing bit) in Eq. ([9]), the 
self-energies involving b must vanish. From this feature 
of the DM field theory, we can derive nonperturbative 
relations among correlation functions as 

JGaait) + G„g(0 + Gacit) - 0, (25) 

JGcait) + GJt) + G,,cit) = (26) 

from the (^',^1'') — (6, a) and (&, c) component of the SD 
equation, respectively. 

In order to study the time evolution of Gaait), we look 
at the (^',^') — ib,a) component, which yields 

Gbait) + TGcait) ~ JAaoGaait) " 2lTTGf,Jt) 

= Fiait)^ (27) 

where 

Fiait) (28) 

= / ds {Sgjt - s)Gaais) + ^l^it - s)G,a(s)} 

^ — oc 

-I f ds {l]gg(i-s)GgJs) + Sjg(i-s)Gga(s)}. 



Using Eqs. ([12]) and (|T4)) . we can rewrite this equation 

for t > as 

Gaait) ~ JTG^ait) + Aa^G aait) = - J F^^it) , (29) 

Now in order to investigate the possible ENE transition 
in this model, we consider the < — > oo limit. Let us 
assume that all the other correlation functions except 
Gaait) vanish in the t ^ oo limit. To the one-loop or- 
der of the perturbation expansion, only the self-energy 
Tiyj^it) can be nonvanishing in the t — > cx) limit due to 
the diagram shown in Fig. [21 since it is proportional to 
(Gaa(O)^- Therefore, the nonvanishing contributions to 
F-^^ioo) come from the first and the third terms on the 
right hand side of Eq. ((28l) . Among the terms on the left 
hand side of Eq. only the third term is nonvanishing 
in this limit. Using Eq. (fT4)) . we therefore have 



JAa^Gaa 



,(oo) = ctG, 



aa(oo) + -Sgg(cX)) 



ds Gads), 






FIG. 2: The one-loop diagram for Egg that may have a non- 
vanishing contribution in the t —> oo hmit. 



where a = ds iE|;^(s) is a finite real number. The 
relation between Gacit) and Gaait) can be obtained from 
the (^f, ^') = (c, a) component of the SD equation, which 
is given by 



- jGaait) + aoGacit) ^ Fcait), 



where 



Fcait) 



(31) 



(32) 



-i / ds {l^cait - s)Gaais) +^ccit ~ s)Gcais)} 

J — OO 

-I f ds {I],^(t - s)Gg„(,s) -f Ege(t - s)GM] ■ 



Now let us suppose that Fcait) can be taken to zero for 
some reason, then, by inserting the expression for Gacit) 
from Eq. ((3T|) into Eq. pO)) . we obtain an equation for 
Gaa (oo), which may have a nonvanishing solution signal- 
ing an ENE transition. This is essentially what happens 
in the field theory with linear FDR as we will see in 
the next subsection. In the DM field theory, however, 
the presence of the first term on the right hand side of 
Eq. ([32]) spoils this scenario. In fact, as i ^ oo, Fcait) 
approaches 



Gc 



/>OC 

i(oo) / ds (-i)Eca(s), 



(30) 



which is nonvanishing by assumption. Then the integral 
in Eq. ([50)1 becomes ill-defined and we are forced to aban- 
don the assumption of the non-zero Gaaioo). This finding 
is consistent with the recent nonpertubative proof by Das 
and Mazenko [21] that the FNH will full spatial depen- 
dence does not support a sharp ENE transition. We note 
that the absence of the ENE transition in our model is 
directly related to the presence of the self-energy Ega in 
our model. This is also similar to the result of Ref . [13] , 
where the self-energy coupling the hatted velocity field 
and the density field plays a crucial role in removing the 
sharp transition. 



B. Field theory with linear FDR 

The SD equation in this case is given similarly to 
Eq. (|24p but now with the component $ taking eight 
different field variables, {(pi, So there are a lot more 
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equations to consider in this approach. The self-energies 
are, however, related to each other through many FDR, 
which can be obtained by applying the time reversal in- 
variance to the SD equation. Here we list only the rele- 
vant FDR among the self-energies to our discussion. For 
(/) = a, 6, 77 or 9, we have 



(33) 



As in the DM approach, we study the time evolution 
of Gaa{t) and its infinite-time limit for the possible ENE 
transition. From the ((f>, $') = (6, a) component of the 
SD equation, we have 

Gbait) + —Gbait) - JAaoGaait) - 2iTTGi,Jt) 

ao 

+TG,,a{t) - JaoGeait) = (34) 



where 



FkJt) ^ -I [ ds ^S^^,(t-s)G^,,(s) 
"'-00 ^, 

-^ f ds J2j:~^^,{t~s)G^,Js). (35) 



Using the FDR in Eqs. ([22]) and ((33l) . we can rewrite the 
above quantity (multiplied by T) as 

TFf,^{t) = [Eg,0(AG,„+G(,a)](i) 

+ [Sgg® ( — Gf,a+G^a)](i) 

ao 

- [Sg-®Gba](i)-[Sg,^®Gaa](t), (36) 

where the convolution between two function f{t) and g(t) 
is defined by 



[f<S>gm= / ds f{t-s)g{s). (37) 
Jo 

As in the DM approach, from the simple form of the 
equation for a{t) and Eq. (j20p . we can derive a nonper- 
turbative relation, namely. 



JGbait) + Gaa{t) = 0. 



(38) 



We now consider the t — > oo limit. If we assume as 
in the previous subsection that all the other correlation 
functions except Gaa{t) vanishes in this limit, then to 
the one-loop order, the only nonvanishing diagram in this 
limit is again the one in Fig. [21 The contribution from 



this diagram to S^g(i) is -2J^ A^Gl^{t). Taking the t 
oo limit in Eq. ([M|) and using Eqs. ([M|) and ([55)) . we 
obtain 



JAaoGaaioo) = AaGaaioo) 

1 



where 



1 /■°° 



(40) 



We therefore have 



Gaa (oo) 



where A = A{1 
parameter by / 



J^ATa^ 



[Gaa{^)-Gaam, (41) 



a/{Jao)). Defining the nonergodicity 

Gaa (oo) /Gaa (0), WC havC 



/ 



1-/ 



where 



C2 



2yl^[Gaa(0)]- 

rial 



(42) 



(43) 



is a dimcnsionless quantity. This is exactly the standard 
MCT equation for the nonergodicity parameter inj;he so- 
called schematic model of the standard MCT [5, 6^ . The 
nonergodic solution / > exists when C2 > 4. 

The origin of the difference between the results of the 
two field theoretic approaches on the existence of an ENE 
transition can be traced back to the terms that are mul- 
tiplied by the self-energy Egg. In both cases, this term is 
expressed as a time integral of a correlation function. In 
the field theory with linear FDR, this correlation function 
is proportional to a total time derivative of Gqq (t) . From 
this, a well-defined equation like Eq. (I42p follows for the 
nonergodicity parameter. In the DM approach, however, 
the correlation function in the integrand is not a total 
time derivative of Gaa{t), but contains an extra contri- 
bution from the self-energy coupling the hatted auxiliary 
field c and the density-like field a. On another level, we 
can understand that the difference comes from the fact 
that, for the one- loop calculation in the field theory with 
linear FDR, only the first order terms are used among 
those in the expression for the fields rj and in Eqs. ([IB)) 
and (|17p . This truncation of the non-polynomial action 
in the perturbation expansion docs not occur in the DM 
field theory. This suggests that the cutoff of a sharp tran- 
sition that the DM approach exhibits at the one loop or- 
der is a kind of nonpertubative information that only an 
infinite resummation in the field theory with linear FDR 
would have an access to. 
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IV. NUMERICAL CALCULATIONS FOR THE 
DAS-MAZENKO FIELD THEORY 

In order to study not just the nonergodicity parame- 
ter, but the full time evolution of Gaa{t), it is desirable 
to have a single equation for the correlation function that 
accounts for the density feedback mechanism. In the orig- 
inal FNH in the DM field theoretic approach [8] , this was 
achieved only in the hydrodynamic limit. It is, however, 
difficult even in the simple toy model to find a time evo- 
lution equation only for Gaa{t)- In this section, we show 
that despite the lack of a complete set of linear FDR, the 
DM field theoretical formulation presents a well-defined 
theory at one-loop order. We do this by constructing 
a closed set of equations for all the correlation functions 
involved and by solving them for Gaa (t) numerically. Be- 
cause of the simple nature of the toy model, especially of 
Eqs. dni), dni), mi), ^ and daHl), we only have five 
independent correlation functions in the DM field the- 
ory. We choose them to be Gaa{t), G^i^{t), Gac{t), G^f^{t) 
and Gcc{t)- AH the other correlation functions can be 
expressed in terms of these five functions. 

If we define the Fourier transform by fito) — 
dt e^^^ f{t) for an arbitrary function /(t), the Fourier 
transforms of these five correlation functions can be writ- 
ten as 



l^^^^|2Ll««MI'(2rr-SggH) (44) 

_ -iJaR{uj) 



-iJTR{uj) 
D{u) ' 

iuj'^ — iJKji{uj) 



where Re denotes the real part and 

Tr{uj) = r + iEg^(a)), 

Kr{uj) = JAaa - 

D{uj) = aR{uj){uj'^ - JKr{uj)) 



(45) 
(46) 
(47) 
(48) 



(49) 
(50) 
(51) 

(52) 



Note that the Fourier transforms of seven self-energies 
are involved in the above equations, which are Eg^(a;), 

S^^(w), Y.ca{ijj), Scc(w), Sgg(w), Sg.(a;) and Ecc(w). At 
the one-loop order of the perturbation theory, these seven 
self-energies are given by functions of the five indepen- 
dent correlation functions thus yielding a closed set of 
equations. It is more convenient to present the one-loop 
self-energies in the time domain. (The detailed derivation 



of the one-loop self-energies in terms of the correlations 
functions will be given elsewhere f28|.) They are given 
by 



~2iaoJ^ATG^~,it)G^j^it) + iaoJTG^^{t)GrAt) 

U^A^Gaa{t)G^l{t) - 2lJATG^l{t)Gac{t), 

1 



-JTG^Q^Sit), 



(53) 



S6,(i) = -lalJ^T (G^s(t))' - ia^JTG^^,{t)GAt) 
+2iaoJ^AT (G^g(t))' - 2JAGaa{t)Gac{t), (54) 

^cait) = -2JAGaait)GAt)+2lJAT{G^^{t)f 

- iTG^^{t)Gcc{t)-iTGac{t)G^i{t), (55) 



^tc{t) = 2iaoJTG^^^{t)G^^^{t)+iTG^i^t)Gac{t) 

- Gaa{t)G,,{t), (56) 



Egg(t) = -2j'A\Gaa{t)f ~2a,J^AT' {G^i{t)y 

+ \alPr'{G^i,{t))\ (57) 



Sg,(i) = ~2iJATGaa{t)G^^{t)+a^JT^G^^,{t)G^^,(t), 



Eg£(t) = -iTGUt)GAt) - (G„g(t))' 



(58) 
(59) 



Note that these expressions are valid only for t > 0. The 
self-energies in Eqs. (|53)) - (j56p vanish for t < due to 
causality. The remaining self-energies satisfy I]gg(— <) = 
Sjg(i), Scs(-t) = Y.,,{t), and I]g-(-t) = -Sg,(t). 

There is one point that requires caution in performing 
a numerical calculation on these types of self-consistent 
equations. We note that Gcc(w) does not decay to zero 
as oj ^ oo. This suggests that there is a delta- function 
singularity in Gcc{t) at short time, which has to be 
treated separately in a numerical calculation. We write 
Gcc{i) = ioi5{t)+ regular terms for some real constant a. 
If we denote by /' and /" the real and imaginary parts 
of a complex function /, respectively, then we can write 



G'U{i^)^a + GT'\-^) 



(60) 



where lim„_oo G"^'°*^^(w) = 0. From Eqs. (US]), (gS]) and 
(15^ . we find that a"^ = ao + lim^^oo S"g(w). At the 
one-loop order, a nonvanishing contribution to E"g(a;) in 
the infinite-oj limit comes from the diagram depicted in 
Fig. [3] We thus have 

lim E^(^) = -a / J^Gaa(^) = -aGaa(O), (61) 
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FIG. 3: The one-loop diagram contributing to E"c(cij) in the 
— > oo Umit 



and 

Gaa(O) 1 1 

2 — 7 ^2 ("2) 

is the initial value of the dimensionless correlation func- 
tion Gaa{t)/aQ. Therefore, the initial value of the cor- 
relation function is determined self-consistently in the 
present field theoretic approach. This is in contrast to 
other approaches [13, 113] where the static limit of the 
correlation functions was used as an input for the ini- 
tial condition. Note that Gaa{0)/o,Q < 1/4 where the 
maximum value occurs when a^a = 2. We believe that 
this peculiar behavior is due to the one-loop perturbation 
theory and that if we consider a higher-loop theory this 
condition will certainly change. 

We now present numerical solutions to the above cou- 
pled equations. We start by evaluating the one-loop self- 
energies using Eqs. (I53p - (l59p from some appropriate ini- 
tial form of the correlation functions (e. g. the bare cor- 
relation functions) given as functions of time. We make 
Fourier transforms of these self-energies and update the 
correlation functions by using Eqs. Then, the 

inverse Fourier transforms are performed to compare the 
input and output correlation functions. This procedure is 
repeated until the convergence is achieved. We find that 
the convergence is achieved in less than 100 iterations in 
most cases. The number TV of mesh points used in the 
time and frequency integrals ranges from 8000 to 26000. 
The cutoffs, Af and A(j, for the time and frequency inte- 
grals, respectively, must be adjusted so that all the five 
correlation functions and the seven self-energies are ac- 
commodated both in the time and frequency spaces. We 
maintain that Af A^^^ ^ Nn to have a consistent numerical 
Fourier transform [29|. As we will explain below, Gaa{t) 
can in general be a relatively long-ranged function, but 
other correlation functions such as G^f^{t) and Gl^g^\t) 
are short-ranged so that we need large A^^ for those func- 
tions. As functions of short and long-ranged are mixed 
in the calculations, both Af and A^^ must be sufficiently 
large. 

For a numerical calculation of the above coupled set 
of equations, we need to put everything in dimension- 
less forms. From Eqs. ([1]) and ([5]), we see that ao/F has 
the dimension of time. Once we put all the correlation 
and response functions in their respective dimensionless 
forms, we find that the self-consistent equations are com- 
pletely described by two dimensionless parameters, k and 



T defined by 

n^^:lp-^\oA, f^l^^^aoT. (63) 

In Figs. [3] and [51 we plot the normalized correlation 
function C{t) = Gaa{t)/Gaa{0) which corresponds to the 
density auto-correlation function in the FNH for various 
values of the parameters k and T. These are compared 
with the corresponding bare correlation functions. From 
these figures, we can see that in general the one-loop 
correlation functions are more stretched in later times 
compared to the bare correlation functions. We perform 
numerical calculations for fixed value of T/k = T/{aQA). 
The analytic expressions for the bare correlation func- 
tions can easily be obtained from the Gaussian part of 
Eq. and we note that the initial value of the bare 

function cia^ {t) is given by G'i°j(0)/a§ = T/k in dimen- 
sionless quantities. As we can see by comparing Figs. S] 
and O the difference between the renormalized correla- 
tion functions and the bare ones is small for small T/k, 
but the nonlinear effects increase with increasing T/k. 
The value of aoa in Eq. ([62)1 is also determined from the 
numerical calculation. For fixed T/k, we find that aoa 
and thus GaaiO)/ciQ is almost constant when we change 
T. As we increase T/k, the initial value Gaa{0)/aQ in- 
creases as well. We find that the numericaljsolutions for 
the self-consistent equations exist only for T/k less than 
some maximum value which is found to be around 0.37 
for the one- loop theory when Gaa{Q)/<il reaches it maxi- 
mum value 1/4. We note that this has nothing to do with 
the ENE singularity which we discussed earlier, since it 
involves the short-time behavior of the correlation func- 
tions. Indeed, as we approach this maximum value, the 
response functions, G^^^it) and G^"^''(i) become increas- 
ingly short-ranged in time. For example, the initial time 
derivative G^g(O^) approaches —oo as T/k approaches 
0.37. We can understand these behaviors by investigat- 
ing carefully the i — *■ 0+ limit of the SD equations ([M]) 
|28| . We believe that the particular values of the pa- 
rameters are specific to the one-loop calculation and will 
change as higher-loop contributions are considered. 

As expected from the discussion in the previous sec- 
tion, C{t) shows a completely ergodic behavior decaying 
to zero as t — > oo for all parameter values in our numer- 
ical calculations. We note that, in the standard MCT 
P, S H, 0] above the ENE transition, the density auto- 
correlation function exhibits a plateau before its eventual 
decay in time. This is not obvious in our numerical re- 
sults. We believe that the present toy model, without 
having a realistic wavenumber dependence, is too sim- 
ple to capture the plateau, if any, in the density auto- 
correlation function calculated from the one-loop approx- 
imation. In order to see if the DM field theoretical ap- 
proach to FNH produces a plateau in a given order of the 
loop expansion, one would have to solve wavenumber- 
dependent versions of Eqs. This remains to be 

seen in the future study. 
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FIG. 4: Normalized correlation function C{t) as a function of 
time t measured in units of ao/F for fixed T/n — 0.1. The 
solid lines are the solutions to the self-consistent one-loop 
equations for f = 0.05,0.025,0.01,0.005,0.0025 and 0.001 
from left to right. The dashed lines are the corresponding 
bare correlation functions. 
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FIG. 5: Normalized correlation function C{t) as a function 
of time t measured in units of ao/F for fixed T/k = 0.2. 
The solid lines are the solutions to the self-consistent one-loop 
equations for f = 0.1,0.05,0.02,0.01 and 0.005 from left to 
right. The dashed lines are the corresponding bare correlation 
functions. 



V. DISCUSSION AND CONCLUSION 

In this paper, we have introduced and studied the toy 
model of FNH of supercooled liquids containing only a 
couple of dynamical field variables a{t) = uq + Sa{t) 
and b{t) without any spatial dependence. We have de- 
veloped two different field theoretic formulations of this 
model first by following the DM prescription and then the 
method involving a set linear FDR. We were also able to 
perform numerical calculations on the coupled equations 
for the correlation and response functions in the DM field 
theory at the one-loop order. 



The major difference between the two field theoretic 
formulations is the way of treating the density nonlinear- 
ities that appears in the problem in the form of l/a{t). 
In the DM field theory this is treated in a simple way 
by introducing a single auxiliary field c(t) and its hatted 
counterpart c(t). On the other hand, in the field theory 
with linear FDR, a couple of auxiliary fields ri{t) and 0{t), 
as well as their hatted partners, are introduced in such a 
way that the dynamical action becomes invariant under 
a set of time-reversal transformations. This manipula- 
tion results in the dynamical action with terms which 
are non-polynomial functions of the main field variables 
Sa{t) and b{t). We note that these terms are proportional 
to ri{t) and 0{t) fields. Therefore, in the renormalized 
perturbation theory at a given order of the loop expan- 
sion, one has to truncate the non-polynomial functions 
at the appropriate order. In some sense, we might say 
that, because of the truncation, the effect of the density 
nonlinearities would not be fully incorporated into the 
field theory at any finite order of the perturbation the- 
ory. We believe that this is one of the reasons why the 
sharp ENE-type transition appears at the one-loop order 
in this formulation. 

This is in contrast to the DM field theory where no 
truncation of the dynamical action is necessary when 
the loop expansion is performed. The loop expansion, 
therefore, has a different meaning from the field theory 
with linear FDR. We have shown that the self-energy 
Sca(i) which corresponds to the one that couples the 
density and the hatted velocity fields plays a key role 
in removing the sharp transition in the one-loop calcula- 
tion. This is consistent with the original finding by DM 
3 and with the recent nonperturbative analysis of the 
FNH [2l|. This fact has often been interpreted as the 
coupling between the current and the density being re- 
sponsible for the ergodicity restoring mechanism. In our 
toy model, the field c(i) is introduced in the DM field 
theory to take care of the density nonlinearities as f)(t) 
and 9(t) are in the other formulation. Even though we 
only perform the one-loop calculation, we can say that 
the self-energy Sca(0 contains the nonperturbative infor- 
mation arising from the density nonlinearities which are 
not present in the field theory with linear FDR. In this 
respect, it may be more appropriate to regard the cutoff 
mechanism resulting from the full nonperturbative treat- 
ment of the density nonlinearities than from the coupling 
between the current and the density. The recent numeri- 
cal calculation '.221 where a direct integration of the gen- 
eralized Langevin equations in FNH are performed also 
demonstrates that the 1/p-nonlinearities (p is the den- 
sity) are playing an essential role restoring the ergodic 
behavior in supercooled liquids. 

It is sometimes discussed in literatures that the DM 
field theory is inconsistent with the FDR. However, a 
set of FDR does hold in the DM field theory, which is 
given in Eq. (fT4|) for our toy model. In the field the- 
ory with linear FDR, a larger set of FDR exists as in 
Eqs. (USD and In the original analysis of 1986 [1], 
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Das and Mazenko used another FDR (linking Gaa and 
Gaa in the toy model notation) in addition to Eq. 
and simplified the equations involved. This is valid only 
in the hydrodynamic limit. In the present paper, we 
do not use such additional simplifications. Instead, we 
keep only Eq. (fT4|) . and take the equations which the 
correlation functions satisfy, Eqs. as a set of 

self-consistent equations, and solve them numerically in 
Sec. IIVI . This program could be generalized to a realis- 
tic situation where the full spatial dependence is present. 
What we have shown in this paper is that the DM field 
theory, viewed as a collection of self-consistent equations 
for the correlation functions at a given order of loop ex- 
pansion, is a well-defined field theoretic approach to glass 
forming liquids. 

There are other field theoretic approaches to the FNH 
than those considered in this paper. A similar toy model 
of the FNH to ours but of different nature has been 
studied by Kim and Kawasaki [23] some time ago. In 
this model, the iV-component of density-like and the 
M-component momentum-like field variables are intro- 
duced without spatial dependence. They consider the 
limit where the numbers N and M approach infinity, and 
find that a sharp transition is present when the condition 
M < N is maintained in the limiting process, while it is 
absent when M = N. In Ref. [qJ, a simpler version of the 
FNH than that of DM was considered, where the sharp 
transition was found to be absent. We note, however, 
that in both cases the effective free energy is quadratic in 
both density and the momentum variables from the out- 
set. There is no need to introduce the auxiliary fields and 
a full set of linear FDR exist in these models. Therefore 
the kind of density nonlinearities discussed in this paper 
is not present and the absence of the sharp transition 
found in these works is probably of a different origin. 

The work by Mayer et al. [30] is another interesting 
zero-dimensional model for glass forming liquids. It is 
in general hard to make a direct connection between the 
projection operator approach, in which Ref. [30] is set. 



and the field theoretical one. It is, however, clear from 
Ref. [13] that a nonperturbative effect that comes from 
considering an infinite number of equations is responsible 
for cutting off the sharp transition. In Ref. [sSj, the sharp 
transition is always present when one considers only a fi- 
nite number of equations. This is related to our finding 
that, when the density nonlinearities are treated in the 
field theory with linear FDR within the loop expansion, 
one has to truncate a non-polynomial function, and the 
sharp transition follows. We might say that the nonper- 
turbative information that cuts off the sharp transition 
is somehow preserved in the DM approach, since one can 
avoid truncating the dynamical action at a given order 
of the loop expansion. 

There are many ways in which the present result can 
be generalized. An obvious generalization is to consider 
a higher-order perturbation theory. Because of the sim- 
ple nature of the model, one can without much difficulty 
construct the higher-loop DM field theory of the model 
and perform the numerical calculations as done in this 
paper. It will be interesting to see how the one-loop re- 
sults, especially the particular initial values of the corre- 
lation functions, get changed when the higher-loop con- 
tributions are considered. We believe that the numerical 
methods developed in Sec. |TV] for treating the coupled 
equations for the correlation and response functions, es- 
pecially those concerning the short-time behavior of func- 
tions, will prove to be useful for an eventual application 
to the full wavenumber dependent FNH in the future. 
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